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is perhaps a natural basis for quantum systems confined to a spherical box in R d , turns out to 
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I. INTRODUCTION 



^ ■ 

a,; 

We contrast two types of confinement for quantum systems, namely confinement in a finite impenetrable box, and 
soft confinement by means of a U-shaped potential. The simplest example is provided by pair of rather different 
problems in dimension d = 1, namely the particle in a box [— L, L] and the harmonic oscillator in R. We use the 

1— i' orthonormal basis {</>i}^ n of the box problem to approximate states of the oscillator. We shall refer to this basis 
as a sine basis since the {4>i\ are scaled shifted versions of the eigenfunctions {\/2 sm(nnx)}'^L 1 for the unit box 
x G [0, 1]. Although the box functions are complete for the Hilbert space L 2 ([— L, L}), they cannot represent the 
oscillator's Hermite functions ipi S L 2 (R) exactly. However, every & is also a member of the Hilbert space L 2 (R). 

— ■ This observation allows us to use the sine functions as variational trial functions for the oscillator. The question 
, remains as to what box size L to use. This is answered by treating L as a variational parameter and minimizing the 

\l ■ upper energy estimates with respect to L. For example, we show in section 2 that by using a sine basis of dimension 



-x 
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N = 50, and optimizing over L, we can estimate the first 10 eigenvalues {1,3,5,..., 19} of the oscillator H = — A 
7— I ■ for d — 1 with error less than 10 -9 . 

In this paper we demonstrate that for problems which are softly confined, or confined to a box whose size is 
greater than or equal to L, the sine functions indeed provide an effective variational basis. In section 2 we study the 
' harmonic oscillator and the quartic anharmonic oscillator in d = 1 dimension. In section 3, we look at spherically 
symmetric attractive potentials in R d , such as the oscillator V(r) — r 2 , the atom V(r) = —1/r, and very singular 
problems V(r) = Ar 2 + Br~ A + Cr -6 , where r g R d : and r = ||r||. Here we employ a sine basis defined on the radial 
interval r £ [a, b], where a and b are both variational parameters. In section 4 we study problems that are themselves 
confined to a finite box p]4Tl| , such as confined oscillators and confined atoms [l|, 0, 0, Q • In section 5 we apply 
the variational analysis two specific many-boson systems bound by attractive central pair potentials in one spatial 
dimension. 



II. PROBLEMS IN R 



In order to work in R, we first consider the solutions to a particle-in-a-box problem confined to the interval [0, 1]. 
By applying the transformation \ = (x — a) / (b — a), we shift the box from the interval [0, 1] to a new interval [a, b]. 



'Electronic address: richard.hall@concordia.ca 
T Electronic address: a'lemusrOlive.concordia.ca 



2 



This gives us new normalized wave functions 

^ Wx)) = /S Sin (^(l^))- (1) 

A special case of this shift is given when the endpoints of the box take the values a = —L and b = L, with L > 0. 
Then we have the explicit wave functions 

Mx(x)) = yjsin (mr (^)) ■ (2) 

We note that the variational basis {^i}^! is a complete orthonormal set for the space L 2 ([a, b]) C L 2 (R), and a 
general element ip S ^ 2 ([ a ; ^]) can be written as the generalized Fourier series 

oo 
i=l 

where Cj = (?/>, 0$) = J ip(x)<fii(x)dx, with i = 1,2,.... This justifies the use of this basis for a variational analysis in 
which the box endpoints {a, b} are to be used as variational parameters. We shall also use a finite basis {<fi n }n=i an d 
include one normalization constraint y] . c„ = 1. This constrained minimization of the expectation value (ipjHip) 
with respect to the coefficients {c„} is equivalent to solving the matrix eigenvalue problem Hv = £v, where 

H = [((f>i, H<pj)] . (3) 

By the Rayleigh-Ritz principle [l2[ for estimating the discrete spectrum of a self-adjoint Schrodinger operator that 
is bounded below, such as H, the eigenvalues £ n of H are known to be one-by-one upper bounds £ n > E n to the 
eigenvalues of H . These bounds can subsequently be further minimized with respect to a and 6, or with respect to L 
in case —a = L = b.. 

Furthermore, to simplify the variational analysis we use the linearity of the operator H, in order to split the matrix 
H in two parts as follows 

H = K + P, (4) 

where K = [{<pi, — A0j)] represents the kinetic energy component, and P = [(4>i,V{x)4>j)] represents the potential 
energy component. The kinetic component will be the same for any potential, in fact, for the basis we have chosen, 
the matrix will be diagonal, where the non-zero elements depend strictly on the variational parameters and have 
analytical exact solution, for example, if we use a box with endpoints {— L,L}, the diagonal elements of K are given 
by i 2 7r 2 /4i 2 , for i = 1, 2, This reduces the total number of calculations required to estimate the eigenvalues of H. 



A. The harmonic oscillator 



It is natural to use the well-known oscillator problem as a test for our variational analysis, since the oscillator is 
not confined explicitly and moreover its eigenfunctions span L 2 (R). We take the scaled one-dimensional harmonic 
oscillator with Schrodinger operator H — — A + x 2 . The solutions to this problem are 

E n =2n+l, (5) 
4> n (x) =c n H n (x)e- x2 / 2 1 (6) 

for n = 0, 1, 2, . . ., where n is the corresponding state, E n the energy of the system, <f> n the wave function, H n the 
Hcrmitc polynomial of order n, and c n is a normalization constant. For this example, we use the basis in Eq.Q, 
with a = —L, b = L, to construct the matrix H. Here, L > is regarded as a variational parameter. We then 
perform a variational analysis using a matrix of dimension N — 50, minimizing the eigenvalues e n over L £ [5.5, 8]. 
We obtain the results shown in Table 1. We note that the absolute approximation error is less than 10~ 9 for the first 
11 eigenvalues. Furthermore, we obtained an error less than 10~ 5 for the first 20 eigenvalues. If we choose a larger 
dimension N for the matrix H, the approximations have a smaller errors, and we can calculate E for higher values of 
n as well, but these results come with a higher computational cost, and can take a long time. 

We can consider the energy levels as functions of the parameter L and fixed N. Figure ([T]) represents the graph of 
the eigenvalues e n of H versus the variational parameter L with fixed N = 50 again, for the first 20 states. We can 
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TABLE I: Approximation of the energy levels of the harmonic oscillator H — —A + x 2 in R. Here, n represents the energy 
state, E the exact solution for the energy, e is the upper bound for E obtained by the variational analysis, with the eigenvalues 
e n of H minimized over the box size, and L is the optimal value obtained. The table shows the energies of the first 12 states, 

71 = 0,1,... ,11 



n 


E 


e 


L 





1 


1.0000000000 


6.86 


1 


3 


3.0000000000 


7.55 


2 


5 


5.0000000000 


7.09 


3 


7 


7.0000000000 


7.14 


4 


9 


9.0000000000 


7.61 


5 


11 


11.0000000000 


7.49 


6 


13 


13.0000000000 


6.85 


7 


15 


15.0000000000 


7.07 


8 


17 


17.0000000000 


7.27 


9 


19 


19.0000000000 


7.43 


10 


21 


21.0000000003 


7.46 


11 


23 


23.0000000017 


7.49 
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FIG. 1: Graph of the eigenvalues e of the matrix H as functions of variational parameter L for fixed N = 50. 

see that these graphs are U shaped and flat near the minima. If N is large enough, the ?7-shaped graphs become even 
flatter, this means that if we take any value of L in this flat region, we will end up with good approximations for the 
energy levels. 

We note that for all calculations in this work, we use the computer algebra software Maple. The advantage of using 
a program such as this is that it does many of the calculations exactly by using symbolic mathematics; it is only at 
the end that it resorts to numerical algorithms to solve, for example, the problem of finding the eigenvalues of the 
Hamiltonian matrix. This minimizes the error obtained in such calculations. 



B. The quartic anharmonic oscillator 

The quartic anharmonic oscillator is another problem in quantum mechanics that has attracted wide interest since 
Heisenberg studied it in 1925. We consider the special case given by the Hamiltonian H = —A + x 2 + x A . Simon 
[l3| wrote an extensive review of this problem and Banerjee et al [lJI calculated the eigenvalues of H using specific 
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scaled basis depending on the harmonic properties of the corresponding eigenfunctions. Analogously to the previous 
section, we have performed a variational analysis, in this case using a matrix of dimension N = 20 and minimizing 
the eigenvalues e„ over L £ [3, 4]. The results are exhibited in Table 2. 



TABLE II: The energy levels for the quartic anharmonic oscillator H = — A + x 2 + x 4 in dimension R. n represents the energy 
state, E represents the accurate energy values obtained in [l4|. e is the upper bound to E obtained using the present variational 
analysis in a basis of size N = 20. The eigenvalues e n of H were minimized over the box size, and L is the optimal value 
obtained. The table shows the first 6 states, n = 0,1, ... ,5. 



11 


E 


e 


L 





1.3923516415 


1.3923516415 


3.4 


1 


4.6488127042 


4.6488127042 


3.4 


2 


8.6550499577 


8.6550499586 


3.4 


3 


13.1568038980 


13.1568038994 


3.7 


4 


18.0575574363 


18.0575574558 


3.4 


5 


23.2974414512 


23.2974415625 


3.4 



We note that with a basis of size only N = 20, the approximation error is of the order of 10 9 for the first five 
states, and then it grows. This problem is solved by taking a larger N in order to reduce the error. 



III. PROBLEMS IN R d 



In order to work in higher dimensions where d > 1 , we need to transform the problem from cartesian coordinates into 
a more suitable system. This approach has been studied in depth by Sommerfeld (HJ. We let x — (x% , . . . , Xd) € R d 
and transform it into spherical coordinates obtaining p — (r,6±, . . . ,9d-i) where r = ||x||. Then the wave function 
will now be given by 

V{p) = 1>{r)Y l (0i,...,6 d -i), 

with ip(r) being the spherically symmetric factor, and Yg the hyperspherical harmonic factor, where £ = 0,1,2,.... 

Given a spherically symmetric potential V(r) in a e?-dimensional space, using the above tools and following, for 
example, the work by Hall, et al. [16|, we get the following radial Schrodinger equation 

d 2 V> d-ldip 1(1 + d-2) . . . _ . 

--JT -J^+ 2 L 1> + V(r)4>=Ety. (7) 

ar z r ar 

Defining the radial wave function 



we rewrite equation Q as 



with effective potential 



R(r) = r^-WVW, R(0) = 



^ + UR = ER, (8) 



g( ,)-y W + <* + '- 1 ff < + '-'> . (9) 



This analysis allows us to work in higher dimensions whenever we consider spherically-symmetric potentials. 



A. The harmonic oscillator 



Using the transformation above, we can work with the harmonic oscillator in higher dimensions, d > 2. A radial 
Schrodinger operator is given now by H = — 4-% + U(r), where U{r) is defined as in Eq. (JSJ and V(r) — r 2 . The 
energy values for this problem are given by 



E n u = 4n + 2£ + d-4, 



(10) 



5 



where i — 0, 1, . . . denotes the angular- momentum quantum number for the d-dimensional problem. The effective 
potential for this problem has a weak singularity and we have found that the variational basis Eq. (fTJ) is suitable 
for such problems, with a = fixed and b > as the remaining variational parameter. However, we do find some 
difficulty in dimension d = 2 when I = : for this specific case, we obtain the effective potential U(r) = r 2 — l/4r 2 . 
The singular term makes the potential tend to — oo when r is close to as shown in Fig. ([2J. This is not an inherent 
feature of the problem but indicates a failure of the effective potential representation when d = 2 and I = : the 
solution to the difficulty is simply to use Eq. (7) as the radial differential equation for this particular case. 



100 



-100 




FIG. 2: Graph of the effective potential U(r) = r 2 - l/4r 2 



We approximate the energy values for the harmonic oscillator in dimensions d = 3, 4, 5 and quantum number 
I = 0,1, 2, 3. We use a matrix of dimension N = 40 and minimize the eigenvalues of H over L £ [3, 12]. The results 
are shown in table . 

If we increase the dimension of the matrix, we see that the error in the calculations decreases, although the 
computer time spent increases considerably. Another example is that of approximating the energy values for the 
harmonic oscillator in dimensions d — 3,4,5 and quantum number I = 0,1,2,3 this time for a larger dimension N. 
We used a matrix of dimension N = 500 and minimized the eigenvalues of H over L £ [4, 12]. The results are shown 
in Table HV1 



B. The hydrogenic atom 

We consider now a special case of the hydrogenic atom in dimension d = 3, that is to say a Schrodinger operator 
given by H = —-^pz + U(r), with U{r) as in Eq. © and V(r) = — — . The energy levels for the model hydrogenic 
atom in this case are given by 



E n £ = 



4(n + ey 



(ii) 



where £ = 0, 1, 2, . . ., and n = 1,2,3,.. .. Since this problem is weakly singular, we use the same basis as in the 
previous example. We calculate approximations to the energy values for the case when e = 1, using a matrix of 
dimension N = 250 minimizing the eigenvalues of H over L £ [3, 190]. The results arc shown in Table IVl 

We see here that the approximation error is larger than 10~ 4 . There are two problems that arise in this analysis. 
First, computations are very slow in this problem due to its singular nature and the number of calculations needed. 
Second, the hydrogen atom has energy levels that are squeezed together as n grows; meanwhile its wave functions are 
very spread-out and quite different from those of the particle-in-a-box problem. This confirms what we would expect 
on general grounds, that the sine basis is not suitable for unconfined atomic problems. 
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TABLE III: Approximation of the energy levels of the harmonic oscillator in dimensions d — 3, 4, 5. n represents the energy 
state, E the exact solution for the energy given by Eq. (|10p . e is the upper bound to E obtained using the present variational 
analysis. The eigenvalues e„ of H are minimized over L. 



d 


£ 


n 


E 


e 


L 


3 





1 


3 


3.00000000 


6.00 






5 


19 


19.00000001 


7.00 






10 


39 


39.00000001 


8.75 




1 


1 


5 


5.00007348 


4.50 






5 


21 


21.00167944 


6.25 






10 


41 


41.00907276 


7.75 




2 


1 


7 


7.00000000 


6.00 






5 


23 


23.00000001 


7.75 






10 


43 


43.00000001 


9.25 




3 


1 


9 


9.00000001 


6.00 






5 


25 


25.00000076 


7.25 






10 


45 


45.00002070 


8.50 


4 





1 


4 


4.00073469 


4.25 






5 


20 


20.00745550 


6.00 






10 


40 


40.02454449 


7.50 




1 


1 


6 


6.00000262 


5.00 






5 


22 


22.00011370 


6.50 






10 


42 


42.00094014 


8.00 




2 


1 


8 


8.00000002 


6.00 






5 


24 


24.00000248 


7.25 






10 


44 


44.00004592 


8.50 




3 


1 


10 


10.00000000 


6.00 






5 


26 


26.00000008 


7.50 






10 


46 


46.00000274 


8.75 


5 





1 


5 


5.00007348 


4.50 






5 


21 


21.00167944 


6.25 






10 


41 


41.00907276 


7.75 




1 


1 


7 


7.00000000 


6.00 






5 


23 


23.00000001 


7.75 






10 


43 


43.00000001 


9.25 




2 


1 


9 


9.00000001 


6.00 






5 


25 


25.00000076 


7.25 






10 


45 


45.00002070 


8.50 




3 


1 


11 


11.00000001 


6.00 






5 


27 


27.00000000 


8.00 






10 


47 


47.00000001 


10.00 



C. Some very singular problems in i? 3 

Problems involving highly-singular potentials are difficult to solve exactly, but they often provide soft confinement 
and may be expected to yield to a variational analysis in the sine basis. Test problems are provided by quasi exactly 
solvable problems. By this is meant that it is possible to find a part of the energy spectrum exactly provided that 
some parameters of the potential satisfy certain conditions. Dong et al. fl7| and Hall et al [l6j studied the potential 

V(r) = Ar 2 + Br- 4 + Cr' 6 

in d = 3. For this work we assume the case where A = 1, C > and £ = 0. Then, for this anharmonic singular 
problem we have the explicit Hamiltonian operator defined by 
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TABLE IV: Approximation of the energy levels of the harmonic oscillator in dimensions d = 3, 4, 5. n represents the energy 
state, E the exact solution for the energy given by Eq. (|1UI) . e is the upper bound to E obtained using the variational analysis. 
The eigenvalues e n of H are minimized over L. This table shows specific examples of energy values for the quantum numbers 
I = 0, 1, 2, 3 and energy states n = 1, 5, 10. Note that the approximation error has diminished compared with those of Table 3. 
In the worst case it is of the of 10~ , while in others cases the upper bounds are almost exact up to numerical accuracy. 



d 


£ 


n 


E 


£ 


L 


3 





1 


3 


3.00000000 


5.4 






5 


19 


19.00000000 


7.4 






10 


39 


39.00000001 


8.7 




1 


1 


5 


5.000000119 


5.4 






5 


21 


21.00000234 


6.9 






10 


41 


41.00001129 


8.4 




2 


1 


7 


7.000000002 


6.0 






5 


23 


23.00000000 


9.0 






10 


43 


43.00000001 


9.0 




3 


1 


9 


9.000000003 


6.0 






5 


25 


25.00000000 


8.1 






10 


45 


45.00000000 


10.0 


4 





1 


4 


4.000009387 


4.7 






5 


20 


20.00008831 


6.5 






10 


40 


40.00027064 


8.1 




1 


1 


6 


6.000000003 


7.7 






5 


22 


22.00000002 


7.4 






10 


42 


42.00000014 


8.8 




2 


1 


8 


8.000000002 


6.0 






5 


24 


24.00000001 


7.4 






10 


44 


44.00000000 


10.0 




3 


1 


10 


10.00000000 


6.6 






5 


26 


26.00000000 


7.6 






10 


46 


46.00000000 


10.0 


5 





1 


5 


5.000000119 


5.4 






5 


21 


21.00000234 


6.9 






10 


41 


41.00001129 


8.4 




1 


1 


7 


7.000000002 


6.0 






5 


23 


43.00000001 


9.0 






10 


43 


43.00000001 


9.25 




2 


1 


9 


9.000000003 


6.0 






5 


25 


25.00000000 


8.1 






10 


45 


45.00000000 


10.0 




3 


1 


11 


1.00000000 


6.5 






5 


27 


27.00000000 


7.6 






10 


47 


46.99999999 


10.3 



r — d o B C „, 
ar r r° 



The exact solution for the ground state is given [l6| by 



E a =4+-J=, (13) 



C 

subject to the constraint (2\fC + B) 2 = C(l + &\/~C). In order to test the sine basis by using a variational analysis 
for this problem, we considered the exact solutions for the ground state energy in two particular cases: first when 
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TABLE V: Approximation of the energy levels of the hydrogen atom in dimension d — 3, for quantum numbers I = 0, 1,2. 
Note that the variational parameter that minimizes the upper bound tends to be very large for all the states. The error is of 
the order of 10 -4 in the 'best' cases. 



1 


n 


E 


e 


b 





1 


-0.2500000000 


-0.2499790730, 


17.5 




2 


-0.06250000000 


-0.06246859682 


40.5 




3 


-0.02777777778 


-0.02773301831 


70.5 




4 


-0.01562500000 


-0.01556528040 


107 


1 


1 


-0.06250000000 


-0.06231120892 


33 




2 


-0.02777777778 


-0.02747649731 


60 




3 


-0.01562500000 


-0.01526320869 


94 




4 


-0.01000000000 


-0.009656788911 


143 


2 


1 


-0.02777777778 


-0.02777640178 


75 




2 


-0.01562500000 


-0.01561644406 


108 




3 


-0.01000000000 


-0.009970374676 


146 




4 


-0.006944444444 


-0.006872824074 


189 



A = B = C = 1 and second, when A = 1, B = C = 9. Since these problems are highly singular, and we are 
considering radial functions, we need to consider two variational parameters, namely the boundaries of the basis 
interval, [a, b]. Thus, we use the basis given by Eq. ([T]) to obtain the matrix H. In this case we need to minimize de 
eigenvalues with respect to a > and b > 0. For A = B = C= lwe have the potential 

V(r) = r 2 +r~ 4 + r~ 6 . 

The ground state energy is given by Eq = 5. We used a matrix of size N = 100 and found that the best result was 
the approximation Eq — 5.00000003, with minimizing parameters a — 0.01 and b = 5.2. For the case where A = 1 and 
B = C = 9 we now have the potential 

V(r) = r 2 + 9r~ 4 + 9r~ 6 

the ground state energy is given by Eq = 7. And our approximation is £o = 7.00000110, where N = 100, and the 
variational parameters that give the minimum value are a = 0.01 and b = 5.1. Even if we have a singular problem, if 
its potential is [/-shaped, we can get upper bounds for the energy levels with a small error. For the sine variational 
basis, the approximations obtained for the upper bound have smaller errors than some of the accurate calculations 
obtained in the references mentioned above. 



IV. CONFINED QUANTUM SYSTEMS 

We can think of this variational approach as if we were confining the system we wish to study in a box, in fact, 
the same box as the particle-in-a-box problem that generates the basis. We need only choose the optimal size to find 
the best approximations to the energy levels. This opens up the possibility to study confined systems themselves, 
provided the basis box size L is less than or equal to the size B of the confining box. The study of these confined 
quantum systems has been of interest in recent years, for example in the early work of Aguilera- Navarro el al [H, 
Michels [f| , Ciftci et al , Al- Jaber 0] , Fernandez and Castro [If. The sine basis yields upper bounds for the energy 
eigenvalues for all L < B. However, we found best results when L = B. This is because the box confinement was 
dominant for the problems studied. Clearly, with potential confinement and a very large B, using an L less than B 
would be advantageous, as it is for unconfined problems. 

A. The confined oscillator 

The confined oscillator was studied by Aguilera-Navarro et al Q who also used the sine basis, with basis box size 
equal to that of the confining box, L = B. We confirm their results, as shown in Table ([VI I) for a box size B = 0.5. 
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TABLE VI: Approximation of the energy levels of the harmonic oscillator H — — 5A+ ^x 2 in dimension d = 1. For the state n, 
E n is a highly-accurate solution for the energy provided by Aguilera-Navarro et al. 3, e n is the upper bound for i? n obtained 
using the present variational analysis with basis size N = 250, and L = 0.5. This table shows the energies of the first 12 states 
n = 0,1,.. .,11 



n 


E 


e 





4.951123323264 


4.951129323244 


1 


19.774534178560 


19.774534179209 


2 


44.452073828864 


44.452073829725 


3 


78.996921150976 


78.996921150748 


4 


123.410710456832 


123.410710456280 


5 


177.693843822080 


177.693843818558 


6 


241.846458758144 


241.846458765623 


7 


315.868612673536 


315.868612686280 


8 


399.760332976128 


399.760332979135 


9 


493.521634054144 


493.521634068796 


10 


597.152524107776 


597.152524136545 


11 


710.653008064512 


710.653008103290 



B. The confined sine-squared potential 

Various confined trigonometric potentials have been studied earlier, for example in Refs. [HI, Ell- We have found 
that these problems can be treated very effectively by a variational analysis in the sine basis. We consider one case 
here, namely the sine-squared potential V(x) given [19| by 




x 



FIG. 3: Graph of the sine-squared potential for Vo = 1- 

By using our variational approach, we immediately obtain the energy eigenvalues exhibited in Table IVIII here, corre- 
sponding to those in Table 1 of For a basis of dimension N — 25, the result differ by at most 10~ 9 . We tabulate 
the relevant results for Vq = 0.1, 1, 5. We have studied both a Hamiltonian Matrix of dimension TV = 10 and another 
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of dimension N = 25: the difference in the results between these two variational bases was found to be of order 10 13 
at most for Vq = 0.1, 1 and the order of 10~ 8 at most for Vq = 5. 

TABLE VII: Approximate energy levels for a confined sine squared potential H = — A + Vbsin 2 (a;) in dimension d = 1. For 
the state n, E denotes the upper bound for the energies obtained using the present variational analysis with basis size N — 25 
and fixed L = ~. The table shows the energies of the first 6 states n = 0,1, ... ,5 



n 


E (V = 0.1) 


E {Vq = 1) 


E (Vo = 5) 





1.024922118883 


1.242428825987 


2.082985293205 


1 


4.049947916808 


4.494793078632 


6.370661125009 


2 


9.050038818610 


9.503664867046 


11.569339156939 


3 


16.050020833189 


16.502081901038 


18.551201398403 


4 


25.050013020839 


25.501302132228 


27.532566336109 


5 


36.050008928573 


36.500892873766 


38.522331587359 



C. The confined atom 

In the case of the unconfined hydrogen atom we found that we needed ever bigger boxes for each following state 
because the wave functions are spread-out. However, the present variational basis is very appropriate for the analysis 
of confined problems themselves. A hydrogen atom confined to a spherical box has been studied by Varshni Q and by 
Ciftci, Hall, and Saad [|[ . In Q , the authors found exact solutions for the confined problem given by the Schrodinger 
equation 

- + ( £ -^±V - -) i>{r) = E^r), (14) 

ar z \ r z r J 

with boundary conditions t/j(0) = t/j(b) = 0, and A > 0. These exact solutions are special for the 3-dimensional case. 
For different quantum numbers, there are specific radii of confinement for which exact solutions are known. These 
problems provide ideal tests for the effectiveness of the sine basis. Details of these exact solutions may be found in 
(9[. We obtain the results shown in Table IVlTll for A = 1 and the radii b required by the available exact solutions. As 

TABLE VIII: Approximation of the energy levels of a confined hydrogenic atom in dimension d = 3. The angular-momentum 
quantum number is I, n is 1 plus the number of nodes of the radial wave function, b is the radius of the box, E is the exact 
value of the energy, E; is the expression in floating point arithmetic, and e is the variational approximation obtained from the 
matrix H of dimension iV = 250. 



i 


n 


b 


E 


Ef 


£ 





1 


4 


-1/16 


-0.06250000000 


-0.0624999668 


1 


1 


12 


-1/36 


-0.02777777778 


-0.0277777498 


2 


1 


24 


-1/64 


-0.01562500000 


-0.0156250000 


3 


1 


40 


-1/100 


-0.01000000000 


-0.0100000000 





1 


3(3 - \/3) 


-1/36 


-0.02777777778 


-0.0277777466 




2 


3(3 + n/3) 


-1/36 


-0.02777777778 


-0.0277775785 


1 


1 


4(5 - VE) 


-1/64 


-0.01562500000 


-0.0156249729 




2 


4(5 + y/E) 


-1/64 


-0.01562500000 


-0.0156248833 


2 


1 


5(7 -V7) 


-1/100 


-0.01000000000 


-0.0100000000 




2 


5(7 + ^7) 


-1/100 


-0.01000000000 


-0.00999999997 


3 


1 


36 


-1/144 


-0.006944444444 


-0.006944444438 




2 


72 


-1/144 


-0.006944444444 


-0.006944444431 



opposed to what we found in the case of unconfined atomic models, it is clear that the sine basis is very well suited 
to the corresponding confined problems. 
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V. THE iV-BODY PROBLEM 



We show in this section that the sine basis can also be effective for the many-body problem. We consider a system 
of ./V identical bosons that are bound by attractive pair potentials V(xi — Xj) in one spatial dimension. In units in 
which h = 1 and m = \ , the Hamiltonian H for this system, with the centre-of-mass kinetic energy removed, may be 
written 



\»=1 \i=l / / l=i<j 



(15) 



By algebraic rearrangement H may be written in the compact form 

JV 



l=i<2 



JV 



+ V(xi - x ) 



(16) 



If ^ is the exact normalized ground state for the system corresponding to the energy E, then boson symmetry allows 
the reduction [20l [2lj to the expectation of a one-body operator whose spectral bottom, in turn, provides an energy 
lower bound El- We have in general, 



E = (*,7J*) = (JV- 1) 



(17) 



where x = (xi — X2). Thus for the harmonic oscillator V(x) — cx 2 , we find immediately that El = (JV — 1)VJV, 
which result coincides in this case with the known [22|, [23[ exact JV-body solution E = c^(JV — 1)\/JV. In order to 
estimate the ground-state energy from above, we employ a single-product trial function $ of the form 



N [2 



cos(nXi/a). 



(18) 



This wave function vanishes outside a box of volume a N in R N . Before we optimize with respect to the box size a, 
we have in general Ejj = ($, H$), where 



E v = (N - 1) 



/7T\ 2 JV 

V ) + ~2 ( ( ^( x i)0( a; 2)l / (a;i - x 2 ),(f>{x 1 )(f>(x2)) 



If we apply Eq. p^|) to the harmonic oscillator F(x) = cx 2 , we find 



E v = {N - l)min 

a>0 



JVca 2 /l 
4 \ 3 



that is to say, 



(7 



c2^(JV - 1)\ZJV, where A = (7r 2 /3 - 2)i w 1.13572. 



(19) 



(20) 



(21) 



Another soluble JV-boson problem is that of the attractive delta potential V(x) = —cS(x). The exact ground-state 

'48 ( 



energy was found by McGuire [24], |25j and is given by the formula E = — -^c 2 JV(JV 2 — 1). Meanwhile the lower and 



upper bounds we obtain respectively from Eqs. (fl"7f and (fT9|) are given by 

E L = -^c 2 (N - 1)JV 2 < E < --^—c 2 (N - 1)N 2 = E v . 



(22) 



The lower bound of course agrees with the exact solution for N = 2. For other numbers of particles, the estimates, 
just as for the corresponding Coulomb one-particle problem, are weaker than those for the tightly bound harmonic 
oscillator. It is also curious that neither bound manages to reproduce the correct JV dependence exhibited by the 
exact formula of McGuire and Mattis. 
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VI. CONCLUSION 



If we compare the harmonic oscillator H = — A + r 2 with the hydrogenic atom H = — A — 1/r in three dimensions 
we see two very different systems from the point of view of stability and the spatial distribution of the respective wave 
functions. The oscillator is tightly bound and hardly exists outside a ball of radius 6, whereas the atom is loosely 
bound and must be considered out to a radius of 50 or more. It is therefore not surprising that the more-confined 
problem, the oscillator, yields to a variational analysis in terms of the sine basis, but the atom does not. The particle 
in a box is the quintessential confined problem. It generates a basis that at first sight might appear inappropriate for 
more general problems. We have shown that it is in fact very effective for problems that are either confined by the 
nature of the potentials involved or are in any case confined by the given boundary conditions. For a systems of N 
identical bosons interacting by attractive pair potentials, the boson permutation symmetry induces behaviour close to 
that of a scaled two-body problem in which the kinetic-energy term is multiplied by (N — 1) and the potential-energy 
term is multiplied by N(N — l)/2. We show that the ground state of this many-body problem can be effectively 
modelled by a product of particle-in-a-box wave functions optimized over the box size L. 
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